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Abstract. There are observational and theoretical indications that both 
the visible (stars) and the dark matter density distributions in elliptical 
galaxies increase significantly up to the galactic center. I present here 
some analytical results obtained with the aid of self-consistent, spheri- 
cally symmetric two component galaxy models. These results suggest the 
possibility that this similar behavior could be a direct consequence of the 
structural and dynamical constraints imposed by the request of positivity 
of the phase-space distribution function (DF) of each density component. 



1. Introduction 

High resolution observations of elliptical galaxies and bulges of spirals show that 
their (spatial) central stellar densities are well described by a power-law profile, 



where < 7 < 2.4 (e.g., Gebhardt et al. 1996; Carollo & Stiavelh 1998). 
Moreover, it is now commonly accepted that a non negligible fraction of the 
total mass in elliptical galaxies is made of a dark component, whose density 
distribution at large radii differs significantly from that of the visible one. The 
radial distribution of the dark matter in the galactic central regions is less well 
known. In the past a commonly used, centrally flat model for the dark matter 
distribution was the so called quasi isothermal (QI) density distribution (/? = 2 
in eq. [6]), but, more recently, high resolution numerical simulations of dark 
matter halo formation showed that also for the dark matter density distribution 



in the central regions, where /? > 1 (e.g., Dubinsky & Carlberg 1991; Navarro, 
Frenk, & White 1997). Thus, there are indications that both the visible and 
the dark matter density distributions increase significantly up to the center of 
galaxies. A natural question arises: is the origin of this qualitatively similar 
behavior independent for the two density distributions, or it is related for some 
reason? 

From a theoretical point of view it is well known that in any physically 
acceptable multicomponent galaxy model the phase-space DF of each density 
component must be non negative. A model satisfying this minimal requirement 
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is called a consistent model. Thus, the question above can be reformulated as 
follows: can the request of positivity of the DF of each density component in 
multicomponent galaxy models tell us something about the relative distribution 
of dark and visible matter in real galaxies? In order to answer this question 
I studied the problem of the construction and investigation of the DFs of two 
component galaxy models with variable amount and distribution of dark matter, 
and variable orbital anisotropy. 



2. Technique 



For a multi component spherical system, where the (radial) orbital anisotropy 
of each component is modeled according to the OM parameterization (Osipkov 
1979; Merritt 1985), the DF of the density component pk is given by: 



fkiQk) = 



1 d f^'^ dQk d'^T 



VStt^ dQk Jo d-^T y/Qk-^T 



Qk{r) = 1 + 



Pk{r), (3) 



' ak, 



where "^rif) = X^k ^ki''') is the total relative potential, and Qfc = ^^^^/^^^^fc- ^ 
and L are respectively the relative energy and the angular momentum modulus 
per unit mass, rak is the anisotropy radius, and fk{Qk) = for Qfc < (e.g., 
Binney & Tremaine 1987). 

Preliminary information on the model consistency can be easily obtained 
using the following results (Ciotti & Pellegrini 1992, CP; Ciotti 1996, C96; Ciotti 
1999, C99). A necessary condition for the non negativity of each fk is: 



dr 



< 0, < r < oo. 
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If the NC is satisfied, strong and weak sufficient conditions for the non negativity 
of each are: 
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> 0, < r < oo. (5) 



The explored two component galaxy models are made of various combi- 
nations of spherically symmetric density distributions, as the fiat-core mass 
distributions 
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the centrally -peaked spatial density associated to the deprojection of the Sersic 
(1968) profile 



I{R) = I{0)ex.p[-b{m){R/Re) 
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(where the explicit expression for b{m) can be found in Ciotti &i Bertin 1999), 
and finally the 7 models (Dehnen 1993) 
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The mass, scale-length, and anisotropy radius of each component are free param- 
eters. For example, in C96 and C99 the family of two component self-consistent 
galaxy models, where one density distribution follows a 71 profile, and the other 
a 72 profile [(71,72) models], is presented. 

3. Results 

The main results can be summarized as follows. 

In CP, by numerical investigation of the inequalities given in eqs. (4)- (5), 
it was shown that it is not possible to add an isotropic QI or Hubble modified 
halo (/3 = 3 in eq. [6]) to a R^/^ galaxy (m = 4 in eq. [7], de Vaucouleurs 
1948), because their DFs run into negative values near the model center. On 
the contrary, the isotropic R^/^ galaxy was found to be consistent for any value 
of the superimposed halo mass and scale-length (i.e., over all the parameter 
space). At variance with the previous case, the isotropic two component models 
where both density distributions are characterized by a centrally flat profile (the 
Hubble modified + QI models), or by a centrally peaked profile (the R^/^+R^/^ 
models), were found to be consistent over all the parameter space. 

In C96 and C99 the analytical DF for both components of OM anisotropic 
(1,1) and (1,0) models is found in terms of elliptic functions. Moreover, the 
method described by eqs. (4)-(5) is applied analytically to the more general 
family of (71, 72) models. It is proved that for 1 < 71 < 3 and 72 < 71 the DF of 
the 71 component in isotropic (71,72) models is nowhere negative, independent 
of the mass and concentration of the 72 component. As a special application 
of this result, it follows that a black hole (BH) of any mass can be consistently 
added at the center of any isotropic member of the family of 7 models when 
1 < 7 < 3, and that both components of isotropic (7, 7) models (with 7 > 1) are 
consistent over all the parameter space. As a consequence, the isotropic 7=1 
component in (1,1) and (1,0) models is consistent. In the anisotropic case, it is 
shown that an analytic estimate of a minimum value of ra for one component 7 
models with a massive BH at their center can be explicitly found. As expected, 
this minimum value decreases for increasing 7. The region of the parameter space 
in which (1,0) models are consistent is successively explored using the derived 
DFs: it is shown that, unlike the 7 = 1 component, the 7 = component 
becomes inconsistent when the 7 = 1 component is sufficiently concentrated, 
even in the isotropic case. The combined effect of the mass concentration and 
orbital anisotropy is also investigated, and an interesting behavior of the DF 
of the anisotropic 7 = component is found and explained: it is analytically 
shown that there exists a small region in the parameter space where a sufficient 
amount of anisotropy can compensate for the inconsistency produced by the 
7 = 1 component concentration on the structurally analogous but isotropic 
case. 

4. Conclusions 

Some of the general trends that emerge when comparing different one and two 
component models with OM anisotropy, as those investigated in CP, Ciotti, 
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Lanzoni, &; Renzini (1996), Ciotti & Lanzoni (1997), C96 and C99, can be 
summarized as follows. 

1) In the isotropic case, when the two density components (i.e., stars and 
dark matter) are similarly distributed in the galactic center, the corresponding 
models are consistent for any choice of the mass ratio and scale-length ratio. 

2) On the contrary, when the two density components are significantly differ- 
ent in the galactic central regions, the mass component with the steeper density 
is "robust" against inconsistency. The mass component with the flatter density 
in the central regions is the most "delicate" and can easily become inconsistent. 

3) For sufficiently small values of the anisotropy radius, OM anisotropy may 
produce a negative DF outside the galaxy center, while the halo concentration 
affects mainly the DF at high (relative) energies (i.e., near the galaxy center). 

4) The trend of the minimum value for the anisotropy radius as a function of 
the halo concentration shows that a more diffuse halo allows for a larger amount 
of anisotropy: in other words, the possibility of sustaining a strong degree of 
anisotropy is weakened by the presence of a very concentrated halo. 

As a consequence of previous points 1 and 2, one could speculate that in the 
presence of a centrally peaked dark matter halo, centrally flat elliptical galaxies 
should be relatively rare, or that a galaxy with a central power law density 
profile cannot have a dark halo too flat in the center. Thus, the qualitatively 
similar behavior of visible and dark matter density profiles in the central regions 
of elliptical galaxies could be a consequence of phase-space constraints. 

Acknowledgments. This work was supported by MURST, contract CoFin98. 
References 

Binney, J., &: Tremaine, S. 1987, Galactic Dynamics, Princeton University Press 

Carollo, CM., & Stiavelh, M. 1998, AJ, 115, 2306 

Ciotti, L. 1996, ApJ, 471, 68 (C96) 

Ciotti, L. 1999, ApJ, 520, 574 (C99) 

Ciotti, L., & Lanzoni, B. 1997, A&A, 321, 724 

Ciotti, L., & Pellegrini, S. 1992, MNRAS, 255, 561 (CP) 

Ciotti, L., Lanzoni, B., & Renzini, A. 1996, MNRAS, 282, 1 

Ciotti, L., Sz Bertin, G. 1999, this conference 

de Vaucouleurs, G. 1948, Ann.d'Ap., 11, 247 

Dehnen, W. 1993, MNRAS, 265, 250 

Dubinski, J., k Carlberg, R.G. 1991, ApJ, 378, 496 

Gebhardt, K., et al. 1996, AJ, 112, 105 

Merritt, D. 1985, AJ, 90, 1027 

Navarro, J.F., Frenk, C.S., k White, S.D.M. 1997, ApJ, 490, 493 
Osipkov, L.P. 1979, Pis'ma Astron.Zh., 5, 77 

Sersic, J.L. 1968, Atlas de galaxias australes. Observatorio Astronomico, Cor- 
doba 



